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Abstract 

Background: Countries aiming for malaria elimination require a detailed understanding of the current intensity of 
malaria transmission within their national borders. National household sample surveys are now being used to 
define infection prevalence but these are less efficient in areas of exceptionally low endemicity. Here we present 
the results of a national malaria indicator survey in the Republic of Djibouti, the first in sub-Saharan Africa to 
combine parasitological and serological markers of malaria, to evaluate the extent of transmission in the country 
and explore the potential for elimination. 

Methods: A national cross-sectional household survey was undertaken from December 2008 to January 2009. A 
finger prick blood sample was taken from randomly selected participants of all ages to examine for parasitaemia 
using rapid diagnostic tests (RDTs) and confirmed using Polymerase Chain Reaction (PCR). Blood spots were also 
collected on filter paper and subsequently used to evaluate the presence of serological markers (combined AMA-1 
and MSP-I19) of Plosmodium folciporum exposure. Multivariate regression analysis was used to determine the risk 
factors for P. folciporum infection and/or exposure. The Getis-Ord G-statistic was used to assess spatial 
heterogeneity of combined infections and serological markers. 

Results: A total of 7151 individuals were tested using RDTs of which only 42 (0.5%) were positive for P. folciporum 
infections and confirmed by PCR. Filter paper blood spots were collected for 5605 individuals. Of these 4769 showed 
concordant optical density results and were retained in subsequent analysis. Overall P. folciporum sero-prevalence 
was 9.9% (517/4769) for all ages; 6.9% (46/649) in children under the age of five years; and 14.2% (76/510) in the 
oldest age group (> 50 years). The combined infection and/or antibody prevalence was 10.5% (550/4769) and varied 
from 8.1% to 14.1% but overall regional differences were not statistically significant (x^ = 33.98, p = 0.3144). 
Increasing age (p < 0.001) and decreasing household wealth status (p < 0.001) were significantly associated with 
increasing combined P. folciporum infection and/or antibody prevalence. Significant P. folciporum hot spots were 
observed in Dikhil region. 

Conclusion: Malaria transmission in the Republic of Djibouti is very low across all regions with evidence of micro- 
epidemiological heterogeneity and limited recent transmission. It would seem that the Republic of Djibouti has a 
biologically feasible set of pre-conditions for elimination, however, the operational feasibility and the potential risks 
to elimination posed by P. vivox and human population movement across the sub-region remain to be properly 
established. 
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Background 

In Africa, some countries in the North have achieved 
malaria eUmination in the last 30 years (Libya in 1973, 
Tunisia in 1979, and Morocco in 2010) and the majority 
of the region is largely regarded as malaria free despite a 
few disputed residual foci in Algeria and Egypt. Elimina- 
tion, however, poses a major challenge in the majority 
of countries in Africa owing to the intrinsically high 
transmission intensity within each country's national 
borders [1] or in the cases of low transmission areas, 
threats posed by neighbouring high transmission coun- 
tries [2], operational constraints posed by fragile health 
systems that may not reach remote foci of infection [2], 
the sophistication of current health information systems 
to identify all cases of clinical and asymptomatic malaria 
[3], and the ability to meet the immediate increases in 
financing needs that would be diverted from other 
health problems [4,5]. 

A recent systematic analysis of available national data 
has defined possible biological threats (including current 
and historical receptivity of transmission and connectiv- 
ity to other countries) and operational threats (including 
political stability, heath system capacity and access to 
populations) to elimination across 99 malaria endemic 
countries [2]. Within Africa the most technically and 
operationally feasible countries, as judged by a combined 
ranking, were Botswana, Swaziland, Sao Tome and 
Principe, Eritrea, and Djibouti. The Republic of Djibouti 
was considered the African country with the most favor- 
able combinations of biological risk and operational 
requirements to reach elimination and its ranking glob- 
ally was equivalent or better than many countries in the 
Americas. However, several recent publications have 
highlighted that despite the fact that 32 countries have 
declared an elimination ambition [6] few have defined 
the true biological feasibility of elimination [7] and fewer 
have defined a combined biological and operational feasi- 
bility assessment [6,8]. There have been limited opportu- 
nities for countries to assess their technical and 
operational strengths and weaknesses for elimination 
unless undertaken as part of very detailed reviews of 
large amounts of empirical data as implemented recently 
in Zanzibar, United Republic of Tanzania [9]. 

National Malaria Indicator Surveys (MIS) have been 
established largely to provide milestone measurements 
for reaching endemic control ambitions in high to mod- 
erate transmission countries, with parasite prevalence as 
the main transmission indicator. In low transmission 
areas this indicator is not efficient [10,11]. Here we pre- 
sent the results of a recent adaptation of a standard MIS 
to include a national serological survey in the Republic 
of Djibouti to define the extent of transmission and 
explore the possibilities of malaria elimination. 



Methods 

National context 

The Republic of Djibouti is a small country of 23,200 Km^ 
located at the southern reaches of the Red Sea on the 
Horn of Africa, bordered by Somaliland region of Somalia 
in the South, Ethiopia in the West and South and Eritrea 
in the North. It was the last French territory on mainland 
Africa to gain independence in 1977 with important stra- 
tegic economic and political value as a port in the Gulf of 
Aden. The significance of the port at Djibouti for the 
Ethiopian economy was recognized at the end of the 19^^^ 
Century when the Emperor declared Djibouti an official 
outlet of Ethiopian commerce that prompted the develop- 
ment of the Addis Abba-Dire Dawa-Djibouti Railway 
(Figure 1). This rail link continues to be an important 
physical connection to Ethiopia. 

The majority of the population are represented by either 
the Issas (ethnically Somalis and connected to clans across 
the sub-region from Somaliland, the Ogaden, and Yemen) 
or the Afars (ethnically connected to people of Ethiopia). 
Other smaller groups exist including people of Arab origin 
and refugees who migrated into the country during the 
1980 s and 1990's from the war-ravaged neighboring coun- 
tries. The population is spatially over-dispersed with 60% 
of the entire country residing in Djibouti Ville and its sub- 
urbs. Five secondary urban settlements exist outside the 
main city at Arta, Obock, Tadjourah, Ali Sabieh and Dikhil 
(Figure 1) and are inhabited by 10% of the population; the 
remaining 30% of the population live in villages or noma- 
dic pastoralists scattered across 90% of the country's semi- 
arid, volcanic land surface in the hinterland. The last 
census was undertaken in 2009 and it was estimated that 
there were 818,159 people [12]. The country is divided 
into six administrative regions. 

There are few rivers but one tributary of the Awash River 
feeds the inland Lake Abbe on the border with Ethiopia 
and seasonal Wadis feed the extremely saline Lake Assal 
which is located west of the inlet water of Ghoubet Kharab, 
and is the lowest point of Africa (Figure 1). Rainfall is 
sparse and erratic with average annual precipitation of 
130 mm with the short rains in March to May and the 
long rains in July to September. The average temperature 
in the country is 30°C with the hot season from May to 
September during which the north eastern monsoon 
blows. Average humidity is above 60% throughout the year 
around Djibouti Ville and much less inland [13]. 

Malaria history and priorities 2006-2010 

The limited annual rainfall and extreme aridity provides 
a natural barrier to stable, endemic transmission of 
malaria across the majority of the Republic of Djibouti. 
At the turn of the last century only two foci of malaria 
transmission were regarded as significant, Ambouli and 
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Figure 1 Map of the Republic of Djibouti showing inland water, roads, railway, administrative regions and major settlements geo- 
coded in Google Earth and obtained from DEPHA [64], showing location of clusters with persons indentified as sero-positive for P. 
falciparum antibodies and/or RDT-PCR confirmed parasite positive (red) and clusters where no person identified as sero-positive or 
parasite positive (green) overlaid on 1 x 1 kilometre resolution surface of aridity as defined by GlobCover [33]. 



Ganaam, on the outskirts of Djibouti Ville [14]. Between 
1910 and the early 1970's almost no locally acquired 
cases of malaria were reported in the Republic of Dji- 
bouti [15,16]. From 1988 epidemics of malaria began to 
appear among refugees from neighboring countries and 
led to onward transmission among local resident com- 
munities [14,15,17]. A study among 144 febrile indivi- 
duals in 1989 in Ambouli area showed no P. falciparum 
infections but exposure of about 6% to antibodies 
against the merozoite stage of the parasite [18]. Ano- 
pheles arabiensis is now accepted as the dominant vec- 
tor of P. falciparum including areas around Djibouti 
Ville as it expands to the wadis, agricultural areas and 
watering holes around the Ambouli region. Some have 
argued that An. arabiensis was introduced to the area 
via the rail link to Ethiopia [19]. Two short malaria 
peaks predominantly occur during the months of May- 
June and November-December [14]. 



The Programme National de Lutte Contre le Paludisme 
(PNLP) was established in 2006 and a National Malaria 
Strategy was launched the same year which recognized a 
possible "re-elimination" end game. The stated ambition 
was to "interrupt the transmission of malaria and elimi- 
nate the disease and its consequences to the socioeconomic 
development of the country as a major public health pro- 
blem in Djibouti' [16]. Despite an anticipated over-disper- 
sion of malaria risks, the aim by 2010 was ''to ensure 
universal coverage of the profitable interventions of preven- 
tion". These were to be achieved by increasing the propor- 
tion of homes owning 2 insecticide treated nets (ITN) to 
>= 90% or 60% of children below the age of 5 years pro- 
tected by an ITN by 2010; introduction of larvivorous fish 
and Abate for larval control; and urban fogging with insec- 
ticides; and case management with first line recommen- 
ded therapy for uncomplicated malaria as artesunate 
and sulphadoxine-pyrimethamine combination (AS+SP) 
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combined with increased diagnostic capacities in all health 
facilities. In addition monitoring systems of health and 
inter-country and cross-border collaboration were to be 
established [16]. These strategies were intended to reduce 
local transmission by half in the whole territory by 2010 
with a view to return to zero the number of actively 
detected cases of malaria and to prevent new transmission 
from existent cases (stated as a pre-elimination objective). 
The Republic of Djibouti was successfully awarded by the 
Global Fund approximately US $ 4 million during round 6 
in 2007 and approved a further award during round 9 in 
2009 amounting to US $6.6 million to mount pre-elimina- 
tion efforts up to 2014 [20,21]. 

The Djibouti 2008-2009 national MIS 
Sampling 

A population probability sample was used to estimate a 
sample of households and clusters for each of the six 
administrative regions reflecting the differences in urban 
and rural settlement within each region. The estimated 
required sample was 4,000 households in 160 clusters 
representing an estimated 24,000 individuals. The sam- 
pling indicator used to power the survey was the pro- 
portion of children under the age of five years sleeping 
under ITN based on the expected change from the Mul- 
tiple Indicator Cluster Survey findings in 2006 [22], pre- 
dicted household size and age composition, between 
cluster design-effect of 1.5, 10% non-participation rate 
and a precision in the estimate of +/- 8.5% derived from 
a balance between optimized precision and costs [23]. 
The survey sample was drawn from a sampling frame of 
450 urban and rural clusters (defined as census enu- 
meration areas with a minimum of 150 households) pro- 
vided by the Direction de la Statistique et des Etudes 
Demographique (DISED) [24]. The sample was then re- 
adjusted to over-sample regions outside of Djibouti Ville 
by reducing the sampled clusters for Djibouti Ville to 
40% of the overall sample rather than the expected 
probability sample of 60% and redistribute the 20% pro- 
portionately across the other less densely populated 
regions of the interior. For each cluster, a household 
sampling interval was selected from the sampling frame 
such that if a cluster had a total of 75 households, every 
third household was surveyed. The sampling interval 
was updated from the village elders during the day of 
survey and any necessary revisions were made. During 
the survey each cluster was positioned using a hand- 
held Global Positioning System (Etrex H, Garmin Inc., 
USA). 

Survey procedures 

A community sensitization exercise was undertaken a 
week before the survey via print, radio and television 
media and local women's groups encouraging support 
and participation during the survey. Training of all 



survey coordinators, supervisors and investigators in 
survey tools and procedures occurred during a two 
week period in November 2008. The survey was under- 
taken over three weeks between December 2008 and 
January 2009 by eight trained teams comprised of a 
team supervisor, interviewer, nurse and laboratory tech- 
nician. Three of the survey teams were assigned to Dji- 
bouti Ville region and one team for each of the other 
regions. The survey instruments used during the house- 
hold investigations were adapted from standard MIS 
tools developed by the Roll Back Malaria Monitoring 
and Evaluation Reference Group [24]. Important varia- 
tions on the standard MIS included ensuring informa- 
tion was captured on all age-groups on ITN use and 
peripheral blood sampling for malaria parasite and anti- 
body detection in a sub-sample of household members 
of all ages. 

At consenting households a de facto census was under- 
taken, including information on age, sex and use of nets/ 
treated nets and their source. Each household head was 
also interviewed about assets ownership and education. A 
random sample of three household members were selected 
from the household census list using a table of random 
numbers and consenting individuals were asked to provide 
a finger prick blood sample. The finger prick sample was 
used to prepare a rapid diagnostic test (Hexagon Pf 
Malaria Rapid HRPII Test (Human GmbH Inc., Germany) 
or ParaHit Rapid Test (ICT Diagnostics, India) to pro- 
vide immediate treatment using nationally recommended 
AS+SP for all positive individuals; a thick and thin blood 
smear, air dried, stained using 5% Giemsa; and three 
5 mm diameter drops of blood on Whatman 3 MM filter 
paper labelled, prepared and stored at 4°C in dry condi- 
tions using procedures described elsewhere [25]. This 
household sub-sample was also interviewed about travel 
histories in the last six months, duration and destination 
and the presence of fever in the last 14 days and on the 
day of survey. 
Laboratory methods 

Serological measures Serologic markers can be particu- 
larly useful in areas of low endemicity, where it may be 
easier to detect relatively long-lasting antibody responses 
than a low prevalence of malaria infections in the human 
population or infected mosquitoes [26-28]. Discs of 
3.5 mm diameter which contained approximately 3 \A of 
blood were punched from each filter paper blood spot 
and reconstituted in low binding 0.5 ml Corning® 
Costar® 96-well flat bottomed plates (Sigma- Aldrich, St. 
Louis, MO, USA) [25]. A 1:100 dilution of whole blood 
(equivalent to a 1:200 dilution of serum) was made by 
adding 300 \i\ of reconstitution buffer (phosphate buf- 
fered saline [PBS], 0.05% Tween-20, 0.1% sodium azide, 
pH 7.2) to each well. The plate was sealed, rocked at 
room temperature overnight, and stored at 4°C until 
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required for the Enzyme Linked Immunoabsorbent 
Assays (ELISA). 

Antibodies were detected simultaneously to P, falciparum 
merozoite surface protein 1-19 (MSP 1 19) and apical mem- 
brane antigen 1 (AMA-1) by ELISA using procedures 
described elsewhere [25]. Antigens were obtained as full- 
length recombinant proteins (AMA-1) or GST fusion pro- 
teins (MSPI19). Immulon 4 96-well plates (Dynatech 
Laboratories, Chantilly, VA, USA) were coated with 50 (il 
of a combination of P. falciparum AMA-1 (3D7) and P. fal- 
ciparum MSPI19 (Kl/Wellcome) at concentrations of 0.5 
(ig/ml in each well and incubated overnight at 4°C. Wells 
were washed 3 times using wash buffer (PBS-Tween 20, 
0.05%). 150 (il/well of blocking buffer (PBS-Tween 20, 1% 
skimmed milk powder) was added to each well and the 
plates were incubated at room temperature for 3 hours. 
Following three washes with wash buffer, 50 (il/well of a 
1:1000 dilution of sera in sample buffer (PBS-Tween 20, 1% 
skimmed milk powder) was added to each well in duplicate 
and incubated overnight at 4°C. Positive controls (hyper- 
immune serum pool) were added in duplicate to each plate 
in a 4-fold serial dilution from 1:50 to 1:51200. Two nega- 
tive controls (European malaria-negative volunteer serum) 
were added to each plate in duplicate at the same concen- 
tration as the test serum. The plates were washed 6 times 
using wash buffer and incubated with 50 (il/well of horse- 
radish peroxidase-conjugated polyclonal rabbit anti-human 
IgG (Dal<o, Roskilde, Denmark) diluted at 1:5000 in sample 
buffer for 3 hours at room temperature. After a further 6 
washes, the plates were developed for 15-20 minutes in the 
dark by the addition of 50 (il/well of SIGMAFA^I^^ o_phe- 
nylenediamine dihydrochloride (OPD) (Sigma- Aldrich, St. 
Louis, MO, USA) substrate solution. The reaction was 
stopped by adding 25 (il/well 2 M H2SO4 and the absor- 
bance measured at 492 nm. 

Parasitology Slides taken during the survey proved to be 
of inadequate quality on reading at the National Public 
Health Laboratory and at the WHO's Regional Centre for 
Excellence in Malaria Microscopy in Oman [Mauolid, 
unpublished observations]. Therefore confirmation of 
RDT positive and 50 randomly selected negative cases 
was undertaken using Polymerase Chain Reaction (PCR) 
protocols for P. falciparum, 

A 1-3 mm^ spot from the filter paper was excised using 
a fresh scalpel blade per sample and put into a 1.5 ml pre- 
labeled tube at room temperature. The QIAamp® DNA 
Mini Kit (QIAGEN Inc, Valencia, CA) was used to extract 
DNA according to the manufacturer's instructions. For 
the RDT positive samples a second extraction was per- 
formed using the ABI6500 semi-automated extraction sys- 
tem. Five positive control samples were produced by 
spotting a parasite blood culture at 40% haematocrit and 
1% parasitaemia. Five negative controls were made from 
non-infected blood in a similar approach. Real-time PCR 



using a P. falciparum specific probe for the the small sub- 
unit of rRNA (18s) gene was carried out using the forward 
and reverse primers 5'-TGCCGACTAGGTGTTGGATG- 
3' and 5' GCCCCAGAACCCAAAGACTT 3' respectively, 
and a TaqMan Minor Groove Binder (MGB) probe 5'- 
FAM-CATCTTTCGAGGTGACTT-MGB-TAMRA 3' 
labeled with 5'FAM (6-carboxyfluorescein) and 3'TAMRA 
(6-carboxytetramethyl-rhodamine) as the reporter and 
quencher, respectively. 2 (il of DNA was added to each 
reaction in a final volume of 20 ml with 1000 nM each of 
the forward and reverse primers, 250 nM of the TaqMan 
probe, and 10 \i\ of TaqMan Universal master mix 
(Applied Biosystems, Life Technologies Corp, Carlsbad, 
CA). The PCR cycling steps were as follows: 50°C for 2 
minutes, 95°C for 10 minutes, and 45 cycles of 95°C for 15 
seconds and 60°C for 1 minute. All positive controls had 
critical threshold (C^) values below 34 cycles. A cutoff of 
38 cycles was used to define positive samples. 
Analysis 

All survey data were entered by four data entry clerks 
using customized data entry packages in MS Access 
(Microsoft Corporation, USA) which included consis- 
tency checks. A wealth index was then calculated from 
household assets indicators using principal component 
analysis [29]. Data were subsequently analyzed using 
STATA version 11 (Statacorp LP, USA). Sampling 
weights were applied due to the disproportionate selec- 
tion of clusters per strata and households within each 
cluster. An additional weight was developed for the para- 
sitology and serology samples where at least 3 individuals 
were sampled per household. The overall weights of the 
samples were equal to product of the inverse of the prob- 
ability of selection of a cluster, household or individual 
and cluster adjusted and weighted estimates were com- 
puted in STATA using the svy command. Visual display 
of all cluster level attributes was undertaken using the 
GPS coordinates for each cluster and displayed using 
ArcGIS 9.3 (ESRI Inc. USA). Inland water bodies were 
assembled from the Global Lakes and Wetlands Database 
(GLWD) developed by the World Wildlife Fund [30] and 
major, minor roads and railway were configured from the 
National Geospatial-Intelligence Agency [31]. The Glob- 
Cover Land Cover product was downloaded from the 
European Space Agency website [32] . This 300 m resolu- 
tion land cover dataset was derived from a time-series of 
Medium Resolution Imaging Spectrometer images 
acquired from December 2004 to June 2006 [33]. The 
GlobCover surface was used to define clusters located in 
extremely arid rural areas that have less than 4% total 
vegetative cover for more than 10 months in a year and 
are therefore least likely to support vector larval survival 
[34,35]. Distances between sampled clusters and aridity, 
main roads and railway and coastline were computed 
using the Spatial Analyst function in ArcGIS 9.3. 
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All duplicate Optical Density (OD) values > 0.2 which 
differed by more than a factor of 1.5 were rejected. To 
define sero-prevalence duplicate OD values were aver- 
aged and converted to titres using the positive control 
curve on each plate to normalize between plates as 
described below. The mean positive control ODs were 
fitted to a titration curve using a sigmoid 3-parameter 
model in Prism version 5 (GraphPad Software Inc, CA, 
USA). Antibody titres were estimated by the equation 
described in [25]. A Gaussian mixtures model was fitted 
to the sample data to determine the sero-negative and 
sero-positive populations underlying the distribution of 
the data [25], each with their own mean and variance. 
This was done using the MCLUST procedure in the R 
programming language [36]. Antibody titres were loga- 
rithmically transformed to render the component distri- 
butions normal. The mean plus 3 standard deviations of 
the sero-negative group was used as the cut-off for sero- 
positivity. A single cut-off was generated as antibodies 
to both AMA-1 and MSP- 1 19 were measured in a single 
well. The sero-conversion rate (k) was estimated by fit- 
ting a reversible catalytic model to the measured sero- 
prevalence by age in years using maximum likelihood 
methods [26,37]. 

PGR positivity was determined using DNA amplifica- 
tion plots of the samples extracted using both the Qia- 
gen and ABI methods which had good correspondence 
between methods. The PGR amplification plots were 
displayed and where curves crossed the critical thresh- 
old for parasite densities and their controls was used to 
determine sample positivity. 

A multivariate regression model was fitted using the 
STATA 11 svy: logistic command allowing for clustering 
and weighting of the data to assess the potential predic- 
tors of P. falciparum exposure and/or prevalence. The 
association of this outcome measure with age, residence, 
use of ITN, fever prevalence within the two weeks of 
survey, distance to main roads and railway, distance to 
the coast line and aridity, travel in the last six months 
and wealth quintile were examined. For each predictor 
the Odds Ratio (OR), 95% confidence interval (GI), and 
P-values were recorded. All predictors were adjusted for 
the effect of the variation between regions. A backwards 
step-wise approach was used to remove the variable 
with the highest non-significant P-value at each stage of 
the regression until a final reduced model was achieved 
with the variables that showed significant association 
with P, falciparum exposure and/or infection prevalence. 

To assess spatial heterogeneity the presence of hot/ 
cold spots of P. falciparum antibody/infection preva- 
lence was examined using the Getis-Ord G-statistic [38] 
and associated Z-scores were computed for each cluster 
in ArcGIS 9.3 (ESRI Inc. USA) using the Spatial Statis- 
tics tool. A high/low value of the G-statistic indicates 



that high/low values of infection and/or exposure are 
clustered within the study area. To determine the signif- 
icance of these statistics, Z-scores were used. When the 
G-statistic is significant, positive/negative Z-scores indi- 
cate clusters of hot/cold spots. 
Ethical approval 

Ethical approval for the study was obtained from the 
Direction de TEpidemiologie et de ITnformation Sani- 
taire (DEIS) of the Republic of Djibouti. All selected 
household heads were informed of the purpose of the 
study and offered opportunities to refuse participation. 
In addition those providing blood samples were 
explained the procedures and informed written consent 
signed for all those agreeing to participate. 

Results 

Defining parasite infection and exposure 

Among the 7717 household members randomly selected 
for parasitological and serological investigation 566 (7.3%) 
refused participation. A total 7151 consenting individuals 
in 156 clusters provided blood samples for testing the pre- 
sence of infection prevalence using RDT with 42 (0.5%) 
positive for P. falciparum infection (Table 1) and con- 
firmed using PGR. The 50 RDT P. falciparum negative 
samples were also confirmed as negative by PGR. Only 
5605 (78%) samples from 129 (83%) clusters were available 
for serological analysis due to interruptions to survey and 
loss of filter papers. Of the successfully sampled persons, 
4769 filter paper samples provided concordant OD data 
for serological measures of AMA-1 and MSP19 P- falci- 
parum antigen exposure (Table 1). P. falciparum sero-pre- 
valence was 9.9% (n = 517) including 46 children aged less 
than five years (5 infants aged less than 6 months). Sero- 
prevalence generally rose with age reaching a peak of 
14.2% in the age group 50 or more years with a significant 
overall trend (x^ = 68.27, p^ 0.0001) (Table 1; Figure 2). 
The seroconversion rate (A,) was estimated by fitting a 
reversible catalytic model to the measured sero-prevalence 
by age [26,37] and was estimated to be 0.0134 new expo- 
sures per person per year (95% GI: 0.010, 0.018) or about 
13 new P. falciparum exposures per 1000 persons per 
year. When the P. falciparum sero-positivity and infection 
prevalence data were combined for the samples that were 
successfully assayed, overall infection + exposure preva- 
lence was 10.5% with a significant increasing trend with 
age {£ = 67.04,^7= 0.0001) (Table 1). 

Analysis of risk factors and spatial heterogeneity of 
parasite infection and/or exposure 

The univariate regression analysis showed that there was 
a statistically significant increase in combined P. falci- 
parum sero-positivity and/or infection prevalence if indi- 
viduals were from the poorest households compared to 
those from the least poor households (Table 2). No 



Table 1 Summary by age of the prevalence of P. falciparum antibody exposures and/or infections, period prevalence of fever two-weeks from the survey day 
and fever on the survey day and use of insecticide treated nets (ITN) the night prior to survey shown as cluster-adjusted and weighted percentages (95%CI) 
[n/N]. 
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For infection prevalence the denominator is all those individuals who gave a finger prick blood sample (n = 7151). For combined infection and antibody prevalence, fever prevalence and ITN use the denominator is 
the total of samples that had concordant optical density values of P. falciparum antibody AMA-1 and MSP-lig antibody exposure (n = 4769). 
^ Seven individuals had missing age data 
^ 14 individuals had missing age data 

^ An individual was regarded as a positive case if he/she showed presence of P. falciparum infection and/or antibodies in the blood sample. Six of the P. falciparum RDT/PCR positive were also positive for antibodies. 
Three individuals who were RDT positive and PCR confirmed were not successfully assayed. Overall, however, infected individuals showed higher titre responses (mean/median = 196.2/98.8) compared to uninfected 
individuals (mean/median titres = 161.3/92.2). 
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significant differences were observed in prevalence of 
infection and/or exposure by urban versus rural; ITN 
use; previous and current fever status; proximity to arid 
areas, the main roads or railway and the coastline; and 
travel history in the six months prior to survey within or 
outside the country (Table 2). The regions of Dikhil 
(14.1%), Tadjourah (14.0%), Ali Sabieh (11.7%) and 
Obock (11.1%) had the highest combined P. falciparum 
infection and exposure prevalence, but the overall differ- 
ences between regions was not significant (Table 2: = 
33.98, p= 0.3144). The multivariate regression model 
showed only age of individual and wealth of household 
were significantly associated with the combined P. falci- 
parum infection and sero-prevalence (Table 3). An 
increase in age group was associated with about 30% 
increase in P. falciparum exposure and/or infection 
prevalence (OR = 1.3; 95% CI = 1.2 - 1.6; p < 0.001) 
(Table 3). An increase in the household wealth quintile 
was associated with a 20% decrease in likelihood of an 
individual being positive for P. falciparum antibodies or 
infections (OR = 0.8; 95% CI = 0.7 - 0.9; p < 0.001) 
(Table 3). The Getis-Ord G-statistic for spatial clustering 



Table 2 Predictors of combined P. falciparum antibody and infection prevalence shown as cluster-adjusted and 


weighted percentages (95%CI) [n/N] 








P. falciparum antibody exposure and infection prevalence 


, p-value 


Rural vs Urban 


12.8 (7.7 - 18.0) [188/1291] vs 9.9 (7.1 - 12.6) [362/3478] 


7.89, p= 0.2906 


Region 






Ali Sabieh 


11.7 (4.1 - 19.3) [57/454] 




Arta 


10.5 (6.2 - 14.9) [25/222] 




Dikhil 


14.1 (7.0 - 21.3) [109/711] 




Djibouti Ville 


8.1 (4.8 - 11.4) [172/2032] 




Obock 


11.3 (5.6 - 17.1) [68/595] 




Tadjourah 


14.0 (6.0 - 21.9) [119/755] 


33.98 p= 0.3144 


ITN use vs not used 


12.2 (6.5 - 17.8) [100/740] vs 10.30 (7.9 - 12.6 [450/4029] 


2.31, p= 0.4456 


Fever in last 14 days Yes vs No 


9.0 (4.1 - 14.0) [26/254] vs 10.6 (8.2 - 13.1) [524/4515] 


0.63, p= 0.5207 


Fever on the day of survey Yes vs No 


15.5 (6.3 - 24.7) [14/91] vs 10.5 (8.1 - 12.9) [536/4678] 


2.35, p= 0.1780 


Cluster located in arid area Yes vs No 


12.0 (8.4 - 15.7) [343/2727] vs 9.7 (5.8 - 11.6) [207/2042] 


13.83, p= 0.1607 


<0.1 km vs > 0.1 km to arid area 


11.5 (8.2 - 14.9) [362/2987] vs 9.0 (5.7 - 12.3) [188/1782] 


7.44, p= 0.3069 


<0.5 km vs >0.5 km to main road/rail 


10.2 (6.7 - 13.7) [274/2456] vs 11.0 (7.6 - 14.3) [276/2313] 


0.74, p= 0.7549 


<7 km vs > 7 km to the coast line 


9.5 (6.7 - 12.4) [276/2466] vs 12.6 (8.3 - 17.1) [274/2303] 


0.01, p= 0.9910 


Travelled anywhere in the last 6 months yes vs no 


10.5 (7.2 - 13.9) [60/403] vs 10.6 (7.0 - 14.1) [490/4366] 


3.03, p= 0.3020 


Most poor vs Least poor 


15.0 (9.0 - 21.0) [171/1063] vs 6.3 (3.4 - 9.2) [64/973] 


41.83, p= 0.0039 




0 20 40 60 80 

Age (years) 



I Predicted Obseived | 

Figure 2 Age specific sero-prevalence plot of P. falciparum 
antibodies using maximum-likelihood fits from reversible 
catalytic equilibrium model [37]. The vertical axis shows the 
proportion of sero-positive individuals in each age group, the 
horizontal axis shows the midpoint age. Red triangles represent 
observed prevalence and blue fitted line is predicted age 
prevalence curve. X = 0.0134 (95% CI, 0.010 - 0.018) and p = 0.0799 
(95% CI, 0.0532 - 0.1200). 
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Table 3 Multivariate regression results of the predictors of combined P. falciparum antibody exposure and/or infection 
prevalence shown as cluster-adjusted and weighted Odds Ratio (95%CI) and P-value of the start full model and final 
reduced model 





Start Full model 




Final reduced model 


Age category^ 


Odds Ratio (95% CI) 

1.4 (1.2 - 1.6) 


P-value 

0.0001 


Odds Ratio (95% CI) P-value 

1.4 (1.2 - 1.6) 0.0001 


Urban vs Rural 


0.9 (0.4 - 1 .7) 


0.660 




ITN use vs not used 


1 .0 (0.6 - 1 .6) 


0.996 




Fever in last 14 days Yes vs No 


0.9 (0.6 - 1 .5) 


0.792 




Cluster located in arid area Yes vs No 


1 .4 (0.8 - 2.4) 


0.214 




<0.5 km vs >0.5 km to main road/rail^ 


0.9 (0.5 - 1 .5) 


0.631 




<7 km vs > 7 km to the coast line*^ 


0.8 (0.5 - 1 .3) 


0.406 




Travelled anywhere in the last 6 months yes vs no*^ 


1.1 (0.7 - 1.8) 


0.361 




Wealth quintile^ 


0.7 (0.6 - 0.8) 


0.0001 


0.8 (0.7 - 0.9) 0.0001 



At each stage of the regression a predictor variable (one with the highest non-significant P-value) was removed from the model in a backwards step-wise 
approach until all variables in the model had a P-value of less than 0.05. 

^ The age grouping shown in Table 1 was used during the regression. 14 individuals had missing age data. 
^ Individuals were classified into the wealth quintile of their resident household. 
^ Median distances were used 



of P. falciparum antibody/infection prevalence was signif- 
icant for 5 clusters with positive Z-scores (hot spots) and 
34 clusters with negative Z-scores (cold spots) while 90 
clusters returned values that suggest non- significant clus- 
tering (Figure 3). All the clusters of high prevalence were 
located in Dikhil region while 33/34 clusters of low pre- 
valence were in Djibouti Ville Region. The analysis of the 
presence of cold/hot spots of parasite prevalence and 
sero-prevalence separately showed no significant cold/ 
hot spots of the former while the pattern for the latter 
was similar to that of combined parasite prevalence and 
sero-prevalence (data not shown). 

Discussion 

We present here the results from the first national MIS in 
sub-Saharan Africa to include combined measures of Pf 
parasitological and serological prevalence among all age 
groups, undertaken in the Republic of Djibouti. P. falci- 
parum parasite prevalence was 0.5% and P, falciparum 
sero-prevalence to MSP-I19 and/or AMA-1 was 9.9% 
(Table 2). Sero-prevalence significantly increased with age 
and the countrywide P, falciparum sero-conversion rate 
{X) was estimated to be approximately 13 new exposures 
per 1000 persons per year. Among the 517 individuals 
who were positive for P. falciparum antibodies, 46 were 
children under the age of five years suggestive of recent 



malaria exposure. Over 76% of these children were from 
Tadjourah, Dikhil and Obock regions (data not shown). 
Overall sero-prevalence in the oldest age group was below 
15%. In 2008 among 2896 patients with malaria-like symp- 
toms attending government run clinics with laboratory 
facilities were tested for presence of infection, 119 (4.1%) 
were reported positive using microscopy and majority of 
these cases (74%) were seen in Djibouti Ville [PNLP, 
unpublished data] . The combined PGR validated P. falci- 
parum RDT positivity and sero-prevalence showed some 
evidence of between region variations with Djibouti Ville 
showing the lowest rates (8.1%). These regional differences 
were overall not statistically significant. This study there- 
fore shows that the risks of P. falciparum malaria are low 
across the RepubUc of Djibouti but majority of clusters 
surveyed, 102/129 (79%), reported at least one individual 
who was positive for P. falciparum antibodies and/or per- 
ipheral blood infection (Figure 3). Although between 
region variations in exposure was not statistically signifi- 
cant, there appeared to be micro-epidemiological hetero- 
geneity in sero-prevalence with 5 clusters located in Dikhil 
region showing relatively higher sero-positivity (Figure 3). 
Residence, aridity, distance to major road/railway and 
coastline, use of ITNs, recent and present fever histories 
and travel were not statistically significant predictors of 
infection and/or sero-prevalence (Table 2). Increasing age 
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Figure 3 The location of clusters that are P. falciparum antibody/infection prevalence hot spots (red dots, n = 5), cold spots (blue 
dots, n = 34) and those that don't show significant clustering (green dots, n = 90) relative to neighbouring clusters. Significance of 
clustering was analysed using the Getis-Ord G-statistic [38]. Positive Z-scores with significant G-statistic indicated hot spots while negative Z- 
scores indicated cold spots. The G-statistic was computed using ArcGIS 9.3 Spotiol Statistics toolbox (ESRI Inc. USA). 33/34 of the cold spots were 
in Djibouti Ville while all hot spots were in Dikhil region. 



and decreasing household wealth were the only two vari- 
ables significantly associated with increasing P. falciparum 
infection and sero-prevalence independently and in 
combination. 

Sero-epidemiological studies carried out in Tanzania on 
AMA-1 and MSP- 1 19 showed that the two antigens pro- 
vided estimates of seroconversion rates that were indepen- 
dently correlated with entomological-based estimates of 
P, falciparum exposure in the same areas [26,39]. AMA-1 
has highly immunogenic properties reflecting recent expo- 
sure [40,41] while MSP-I19 antibodies saturate over longer 
periods of exposure representing cumulative risks [26]. 
Combining exposure to both antigens in a single ELISA 
test is likely to provide a more robust measure of parasite 
exposure in areas of low parasite prevalence. However 
interpretation of sero-positivity from optical density values 
of titres measured through ELISA has an unknown sen- 
sitivity and specificity and will depend largely on the 



approaches taken to classify positivity. Here we elected not 
to use traditional three standard deviations above estab- 
lished control sera but have used a mixture model method 
which has previously been used for defining sero-positivity 
to measles and rubella [42-44]. By using the population 
itself to define sero-positivity we remove any bias due to 
population-specific immunoglobin levels, nutritional or 
co-infection history. However, defining sero-positivity still 
requires further work to either define their sensitivity, 
eliminate the need for a cut off or improve approaches to 
defining it. 

From a sampling perspective the largest limitation of the 
MIS was that it was powered to detect changes in ITN 
coverage and not parasite exposure. This represents a per- 
ennial problem for national MIS across Africa designed to 
detect increases in coverage measures and not declines in 
exposure. There is a recommendation however by partners 
to measure changing parasite prevalence through national 
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surveys every 3-5 years [45]. Using the period prevalence 
of parasitaemia and/or sero-positivity and the between 
cluster variability from the survey undertaken in the 
Republic of Djibouti in 2008-09 we estimate that one 
would require testing a total population of over 22,000 to 
achieve a +/- precision of 5% within a point estimate. As 
prevalence becomes rare the costs of sampling to define 
parasite exposure therefore increase substantially. 

The field investigations of parasite prevalence included 
two types of HRP-2 RDTs specific only to P. falciparum 
and these, and a random sample of 50 RDT negatives, 
were validated using PGR resulting in exactly similar 
results. It is widely held in the Republic of Djibouti that P, 
vivax constitutes less than 5% of all malaria infections [16], 
however this remains to be investigated properly as vivax 
infections are often missed during routine clinical micro- 
scopy [46] . We have not established P, vivax prevalence 
during the national survey and due to funding constraints 
we were not able to undertake large-scale PCR-based 
investigations of vivax and other rare species on the filter 
paper samples. All neighboring countries connected to the 
Republic of Djibouti also have P, vivax transmission. P, 
vivax is harder to detect, treat and eliminate [47]. The 
national guidelines on malaria case-management for the 
Republic of Djibouti are not clear on recommendations for 
the diagnosis and treatment of vivax malaria [48]. 

In 2009 the Republic of Djibouti submitted a proposal to 
the GFATM for five years support to achieve a malaria 
pre-elimination status by 2015. About 6.6 million USD 
was approved in June 2009 but the agreement is unlikely 
to be signed before mid-2011 (personal communication. 
The Global Fund). Our findings support the biological 
possibility for elimination reported elsewhere [2] but do 
not examine the operational feasibility of elimination. 
Interestingly the GFATM support from round 6 and 
round 9 continues to promote universal coverage of ITN 
[20,21]. ITN coverage was 1.4% among children in 2006 
[23] and rose in 2008-09 to 19.6% among children aged 
less than 5 years, with 41.5% of households owning at least 
one ITN and 14.4% of the entire population using them 
the night before survey [22]. Increasing ITN coverage to 
over 80% would continue to be a huge resource outlay for 
little obvious biological benefit where malaria risk is intrin- 
sically low and spatially over-dispersed. A more logical 
approach might involve targeted vector control with com- 
binations of larval and adult mosquito target interventions 
in response to improved vector breeding and/or clinical 
case detection and cartography [49]. Mathematical models 
suggest that targeting the minority who contribute the 
majority to transmission in low transmission settings is 
likely to yield larger returns [50-52]. Although there was a 
generally limited clustering of antibody/infection preva- 
lence in Djibouti, there are 4 regions (Tadjourah, Dikhil, 
Obock and Ali Sabieh) which contribute to the majority of 



cases. However a migration from blanket coverage to tar- 
geted intervention is predicated on a robust and timely 
surveillance system transitioning from routine clinical 
treatment with incomplete monthly and quarterly report- 
ing to a model based on daily or weekly aggressive, infec- 
tious disease surveillance supported by public awareness 
campaigns and universal parasitological testing. Malaria is 
not a notifiable disease in Djibouti and case reporting is 
embedded within the routine health information system 
which suffers from perennial under-reporting and data 
quality problems typical of most health information sys- 
tems in Africa [53-55]. Surveillance precedents already 
exist in the Republic of Djibouti for polio and cholera and 
new infectious disease surveillance initiatives are being 
developed to track West Nile and Dengue viral infections 
[Amar Abdo, personal communication]. Malaria surveil- 
lance also needs to become an up-regulated disease 
reporting system with the ability to detect, intervene and 
eliminate residual foci of malaria infections as and when 
they emerge. The round 9 proposal includes 0.48 million 
USD to strengthen surveillance but deciding on how best 
to develop and implement malaria surveillance systems for 
elimination has a much weaker evidence-based platform 
than ITN delivery [3]. 

Although travel histories in the last six months were not 
associated with P, falciparum infection or sero-prevalence 
in this study, human population movement is likely to be a 
significant challenge to elimination in the Republic of Dji- 
bouti. Several authors have highlighted the role played by 
the rail route to Ethiopia in the recent re-emergence of 
malaria transmission in Djibouti [14,56,57]. Between 
700,000 and 1.4 million passengers use this rail route each 
year [58]. Residents of Djibouti Ville often migrate during 
the months of May- September to avoid the excessive heat 
to Hargeisa in Somaliland or Dire-Dawa in Ethiopia. The 
road from Djibouti Ville to Addis via Galafi has recently 
been upgraded with funding from the World Bank and the 
European Union to increase transportation of goods and 
people between countries [58]. Djiboutian armed and 
police forces train in Ethiopia and are stationed as part of 
UN peace keeping in Cote Dlvoire. Thousands of foreign 
military personnel are stationed in Djibouti Ville and 
engage in regional travel. Therefore the significance of 
human population movement for the risks of imported 
infectious diseases, including malaria, should not be under- 
estimated and may be an important contributor to disease 
hot spots in the country. Although recent reports show sig- 
nificant reductions in malaria infection prevalence in Ethio- 
pia [59], Eritrea [60], Somalia [61] and Yemen [62], the 
connectivity between the Republic of Djibouti and these 
countries demonstrates the importance of constantly 
reviewing the origin and destination risks against regional 
human population movement to define and contain impor- 
tation risks [63]. Use of mobile phone data to track 
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population movements internally and externally [63]; 
screening arriving passengers for malaria at the border 
entry points and recording of their detailed travel history; 
and redesigning household survey sampling and question- 
naires to capture detailed travel histories of individuals 
tested for malaria infections are possible approaches to 
assessing the impact of human population movement on 
malaria in Djibouti. Work is ongoing to assess the theoreti- 
cal impact of human population movement on malaria 
transmission in Djibouti using a combination of travel his- 
tory data from surveys; air passenger origin and destination 
data; and estimates of malaria risk of origin and destination 
[Deepa Pindolia, Personal Communication]. 

Conclusion 

This study demonstrates the potential of adapting stan- 
dard national MIS to collect information not only on 
the prevalence of malaria infections but also serological 
markers of exposure which are important in determin- 
ing transmission in areas of low endemicity [26]. The 
financial cost associated with the materials required to 
collect dry blood spots for serology and subsequent ana- 
lysis were relatively small. The additional survey time as 
a result of collecting the samples on the filter papers 
appeared minimal. However, it was not possible to con- 
duct a full cost-effectiveness analysis of this adaptation 
of a standard MIS as this was beyond the scope of this 
study. The biggest challenge, perhaps, is less to do with 
actual cost of including a serological component in stan- 
dard MIS but the expertise required to implement 
ELISA to investigate exposure. This serological method 
for malaria is highly specialised and currently implemen- 
ted largely by research groups. Potentially, however, 
training of technicians at national public health labora- 
tories in malaria ELISA, as was done in Djibouti for this 
study, could mitigate this constraint. 

The results of the survey support the apparent low 
case reporting of P. falciparum in the Republic of Dji- 
bouti and at face value the country may appear to be 
poised for elimination, a state enjoyed across the coun- 
try between 1910 and the 1970's [15,16]. P. falciparum 
elimination would be possible if surveillance and 
response to foci of infections were substantially 
improved drawing on experience and systems estab- 
lished for polio and other rare, notifiable infectious dis- 
eases. Detailed investigations into the extent of P. vivax 
transmission are also needed for a comprehensive 
understanding of the possibility of elimination [47]. This 
would require careful planning through a feasibility ana- 
lysis as undertaken recently in Zanzibar [9] and long- 
term uninterrupted political support and financing [8]. 
At present the national malaria control programme has 
managed to secure some funding to support salaries 
from the GFATM round 6 funding for a short period 



but no funds since 2008 to implement any activities 
[Hawa Guessod, Personal Communication]. The signing 
and releasing of funds from round 9 has yet to materia- 
lize although it is anticipated that a larger national sur- 
vey will be undertaken in 2012 during which an RDT 
that can detect both P, falciparum and P, vivax will be 
used in addition to collecting blood samples on filter 
papers for multi-species PCR and ELISA [Hawa Gues- 
sod, Personal Communication]. This and the proposed 
future studies will substantially improve our understand- 
ing of malaria transmission in Djibouti, but without a 
carefully defined plan, staff or funding elimination of 
malaria in Djibouti remains only a biological possibility. 
The threshold of parasite prevalence seen as the bench- 
mark for deciding whether to move on to malaria elimi- 
nation or sustain low endemic control is seen as 1% 
parasite prevalence [3,11]. But to detect this level of pre- 
valence requires large and expensive survey samples. 
The threshold for serological markers remains unclear 
although the absence of exposure among children under 
the age five years would be a good indicator that a 
country is approaching pre-elimination [3]. An alterna- 
tive equivalent index is case incidence of 1 per 1000 
persons at risk [3,26]. To reliably estimate this index, 
there is need for accurate passive and active case detec- 
tion data over several years [11]. This requires a prop- 
erly functioning national malaria surveillance system 
based on a quality assured diagnostic capacity to provide 
the accurate information on whether Djibouti has 
achieved the biological threshold for elimination. 

Acknowledgements 

The Malaria Indicator Survey was implemented by the Programme National 
de Lutte Contre le Paludisme of the Ministry of Health of the Republic of 
Djibouti with funding support from the Global Fund and the United Nations 
Children's Fund. Technical assistance was provided by the World Health 
Organization's Eastern Mediterranean Regional Office of (WHO-EMRO) and 
the Wellcome Trust, UK. We also wish to acknowledge Bart Faber and Ed 
Remarque of the Biomedical Primate Research Centre, Rijswijk, The 
Netherlands, for providing the antigen used in this study. AMN is supported 
by the Wellcome Trust as a Research Training Fellow (#081829). CKM is 
supported by a Wellcome Trust grant to the Kenya Major Overseas 
Programme (#084538). MJM is supported by a Wellcome Trust grant 
(#088634). RWS is supported by the Wellcome Trust as Principal Research 
Fellow (#079080). CJD acknowledges the support of the Wellcome Trust 
(#078925) which also supports JC. The authors are grateful to the support 
and diligence of the 32 field staff notably the supervisors Abdallah 
Mohamed, Abdillahi Ayeh, Mohamed Abakari, Abdoulkader Garad, Mohamed 
Moussa, Idriss Farah and Mahad Ibrahim and to Viola Kirui and Deepa 
Pindolia for help with data cleaning, preparation and map production. We 
thank Dr Tharcisse Barihuta for his advice and facilitation of the survey. The 
authors would like to thank Mohamed Abakari, Ubah Moussa, Abdinour 
Sudan for their help in preparing the ELISA samples at the Laboratoire 
National Sante Publique. The authors are also grateful for comments on 
earlier drafts of the manuscript from Dr Emelda Okiro. 

Author details 

^Malaria Public Health & Epidemiology Group, Centre for Geographic 
Medicine Research - Coast, Kenya Medical Research Institute/Wellcome Trust 
Research Programme, P.O. Box 43640, 00100 GPO, Nairobi, Kenya. ^Centre for 



Noor et al. BMC Infectious Diseases 201 1, 1 1 :1 21 
http://www.biomedcentral.eom/1471-2334/1 1/121 



Page 13 of 14 



Tropical Medicine, Nuffield Department of Clinical Medicine, University of 
Oxford, CCVTM, Oxford 0X3 7LJ, UK. ^Direction de I'Epidemiologie et de 
rinformation Sanitaire, Ministere de la Sante, P.O. Box 438, Djibouti, 
Republique de Djibouti. ^Programme National de Lutte Contre le Paludisme, 
Ministere de la Sante, P.O. Box 1157, Djibouti, Republique de Djibouti. 
^Centre for Geographic Medicine Research - Coast, Kenya Medical Research 
Institute/Wellcome Trust Research Programme, P.O. Box 230, Kilifi, Kenya. 
^Department of Immunity and Infection, Faculty of Infectious and Tropical 
Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, 
London, WCIE 7HT, UK. 

Authors' contributions 

AMN was responsible for study design, data cleaning, analysis, interpretation, 
drafting and production of the final manuscript. MBM was responsible for 
the serological analysis of samples in the laboratory and data assembly and 
contributed to the final manuscript. CKM was responsible for the analysis of 
the optical density data, interpretation of results and contributed for the 
final manuscript. MAO, HHG, lAA contributed to the survey design, data 
assembly and cleaning and contributed to final manuscript. MN was 
responsible for the polymerase chain reaction analysis of all positive and 
selected negative samples and contributed to the final manuscript. JC and 
CJD were responsible for the development of ELISA methods, interpretation 
of results and contributed to final manuscript. MJM contributed to the 
statistical analysis of the optical densities, was responsible for the PCR of 
positive and negative samples and contributed to the final manuscript. RWS 
was responsible for overall scientific management, analysis, interpretation 
and preparation of the final manuscript. All authors read and approved the 
final manuscript. 

Competing interests 

The authors declare that they have no competing interests. 

Received: 14 January 201 1 Accepted: 1 1 May 201 1 
Published: 11 May 2011 

References 

1. Snow RW, Marsh K: Malaria in Africa: progress and prospects in the 
decade since the Abuja Declaration. Lancet 2010, 376:137-139. 

2. Tatem AJ, Smith DL, Gething PW, Kabaria CK, Snow RW, Hay SI: Ranking 
elimination feasibility in malaria-endemic countries. Lancet 2010, 
376:1579-1591. 

3. Moonen B, Cohen JM, Snow RW, Slutsker L, Drakeley C, Smith DL, 
Abeyasinghe R, Rodriguez Lopez MH, Maharaj R, Tanner M, Targett G: 
Operational strategies to achieve and maintain malaria elimination. 
Lancet 2010, 376:1592-1603. 

4. Sabot 0, Cohen JM, Hsiang MS, Kahn JG, Basu S, Tang L, Zheng B, Gao Q, 
Zou L, Tatarsky A, Aboobakar S, Usas J, Barrett S, Cohen JL, Jamison DT, 
Feachem RGA: Costs and financial feasibility of malaria elimination. 
Lancet 2010, 376:1604-1615. 

5. Snow RW, Okiro EA, Gething PW, Atun R, Hay SI: Equity and adequacy of 
international donor assistance for global malaria control: an analysis of 
populations at risk and external funding commitments. Lancet 2010, 
376:1409-1416. 

6. Feachem RGA, Phillips AA, Hwang J, Cotter C, Weilgosz B, Greenwood BM, 
Sabot 0, Rodriguez Lopez MH, Abeyasinghe R, Ghebreyesus TA, Snow RW: 
Shrinking the malaria map: progress and prospects. Lancet 2010, 
376:1566-1578. 

7. Cohen JM, Moonen B, Snow RW, Smith DL: How absolute is zero? An 
evaluation of the historical and current definitions of malaria 
elimination. Malaria Journal 2010, 9:e213. 

8. Feachem RGA, Phillips AA, Targett GA, Snow RW: Call to action: priorities 
for malaria elimination. Lancet 2010, 376:1517-1521. 

9. Zanzibar Malaria Control Program: Malaria Elimination on Zanzibar: A 
Feasibility Assessment. 2009 [http://www.malariaeliminationgroup.org/ 
malaria-elimination-zanzibar-feasibility-assessment], Accessed December 
2010. 

10. Yekutiel P: Problems of epidemiology in malaria eradication. Bulletin 
World Health Organ 22:669-683. 

11. Hay SI, Smith DL, Snow RW: Measuring malaria endemicity from 
intense to interrupted transmission. Lancet Infectious Diseases 2008, 
8:369-378. 



12. Direction de la statistique et des etudes demographiques: Resultats 
preliminaires du deuxieme recensement general de la population et 
I'habitat 2009. Ministere de I'economie des finances et de la planification, 
charge de la privatisation 2010. 

13. Climatetemp. [http://www.climatetemp.info/djibouti/]. Accessed January 
2011. 

14 Rodier GR, Parra JP, Kamil M, Chakib SO, Cope SE: Recurrence and 

emergence of infectious diseases in Djibouti City. Bulletin of the World 
Liealth Organization 1995, 73:755-759. 

15. Manguin S, Carnevale P, Mouchet J, Coosemans M, Julvez J, Richard- 
Lenoble D, Sircoulon J: Biodiversity of malaria in the world. John Libbey 
Eurotext, Paris , English 2008, 140-141. 

16. Programme National de Lutte Contre le Paludisme: Plan strategique la 
lutte contre le paludisme 2006-2010. Ministere de la Sante, Republique de 
Djibouti 2006. 

17. Louis JP, Albert JP: Le paludisme en Republique de Djibouti strategie de 
controle par la lutte antilarvaire biologique poissons larvivores 
autochtones ('aphanius disparet toxines bacteriennes. Medecine Tropical 

1988, 48:127-131. 

18. Fox E, Abbatte EA, Leef M, Mikhail E, Said-Salah Y, Hassan A: Malaria in the 
Djibouti Republic-Results of a serologic survey in Ambouli. Medecine 
Tropicale 1989, 49:159-160. 

19. Rogier C, Pradines B, Bogreau H, Koeck JL, Kamil MA, Mercereau-Puijalon 0: 
Malaria epidemic and drug resistance, Djibouti. Emerging Infectious 
Diseases 2005, 11:317-321. 

20. Global Fund to Fight AIDS, IB and Malaria Round 6 Djibouti proposal. 
[http://www.theglobalfund.org/grantdocuments/6DJBM_1307_0_full.pdfl, 
Accessed January 201 1. 

21. Global Fund to Fight AIDS, IB and Malaria Round 9 Djibouti proposal. 
[[ittp://www.t[ieglobalfund.org/grantdocuments/9DJBM_1828_0_full.pdfl, 
Accessed January 201 1. 

22. Enquete Djiboutienne a Indicateurs Multiple 2006: [[ittp://www.c[iildinfo. 
org/mics3.[itml]. Accessed December 2010. 

23. RBM-MERG: Monitoring and Evaluation Resource Group. Core Household 
Questionnaire. [[ittp://www.rbm.w[io.int/partners[iip/wg/wg_monitoring/ 
docs/mis2005/ccl.pdf], Accessed December 2010. 

24 Programme National de Lutte Contre le Paludisme: Djibouti National 
Malaria Indicator Survey 2008-2009. World Health Organization, Eastern 
Mediterranean Regional Office, Cairo, Egypt 2009. 

25. Corran PH, Cook J, Lync[i C, Leendertse H, Alp[iaxard M, Griffin J, Cox J, 
Abeku T, Bousema T, Ghani AC, Drakeley C, Riley E: Dried blood spots as a 
source of anti-malarial antibodies for epidemiological studies. Malaria 
Journal 2008, 7:el95. 

26. Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SL, Carneiro I, 
Malima R, Lusingu J, Manjurano A, Nkya WM, Lemnge MM, Cox J, 
Reyburn H, Riley EM: Estimating medium- and long-term trends in 
malaria transmission by using serological markers of malaria exposure. 
Proceedings of National Academy of Science, USA 2005, 102:5108-5113. 

27. Okell LC, GUani AC, Lyons E, Drakeley CJ: Submicroscopic infection in 
Plasmodium foldporum-en6em\c populations: a systematic review and 
meta-analysis. Journal of Infectious Diseases 2009, 200:1509-1517. 

28. Bousema T Youssef RM, Cook J, Cox J, Alegana VA, Amran J, Noor AM, 
Snow RW, Drakeley C: Serological markers can detect heterogeneity in 
malaria exposure in an area of very low transmission intensity in 
Somaliland. Emerging Infectious Diseases 2010, 16:392-399. 

29. Filmer D, Pritc[iett L: Estimating wealth effects without expenditure data- 
or tears: and application to educational enrolments in states of India. 
Demography 2001 , 38:1 1 5-1 32. 

30. iehne( B, Doll P: Development and validation of a global database of 
lakes, reservoirs and wetlands. Journal of Hydrology 2004, 296:1-22 
[[ittp://www.worldwildlife.org/science/data/iteml877.[itml], Accessed 
December 2010. 

31. National Geospatial-lntelligence Agency (NGA): [[ittp://eart[i-i nfo.nga.mil/ 
publications/vmapO.[itml], Accessed December 2010. 

32. European Space Agency: [[ittp://ionial.esrin.esa.int/]. Accessed December 
2010. 

33. Arino 0, Gross D, Ranera F, Leroy M, Bic[ieron P, Brockman C, Defourny P, 
Vancutsem C, Ac[iard F, Durieux L, Bourg L, Lat[iam J, Di Gregorio A, Witt R, 
Herold M, Sambale J, Plummer S, Weber JL, Goryl P, Houg[iton N: 
GlobCover: ESA service for global land cover from MERIS. IEEE 
International Geoscience and Remote Sensing Symposium 2007, 2412-2415. 



Noor et al. BMC Infectious Diseases 201 1 , 1 1 :1 21 Page 1 4 of 1 4 

http://www.biomedcentral.eom/1471-2334/1 1/121 



34. 



35. 



36. 



37. 



38. 



39. 



40. 



42. 



43. 



44. 



45. 



46. 



47. 



49. 



50. 



52. 



53. 



54. 



55. 



Omer SM, Cloudsley-Thompson J: Survival of female Anopheles gambiae 
Giles through a 9-month dry season in Sudan, Bulletin of World Health 
Organization 1970, 42:319-330. 

Guerra CA, Howes RE, Patil AP, Gething PW, Van Boeckel TP, Temperley WH, 

Kabaria CW, Tatem AJ, Manh BH, Elyazar IRF, Baird KJ, Snow RW, Hay SI: The 

international limits and population at risk of Plasmodium vivax 

transmission in 2009. PLoS Neglected Diseases 2010, 4:e774. 

Fraley C, Raftery A: MCLUST: Model-based clustering and normal mixture 

modeling including Bayesian regularization. [http://www.stat.washington. 

edu/feseafch/repofts/2006/tf504.pdf]. Accessed December 2010. 

Corran P, Coleman P, Riley E, Drakeley E: Serology: a robust indicator of 

malaria transmission intensity? Trends in Parasitology 2007, 23:575-582. 

Getis A, Ord JK: The Analysis of Spatial Association by Use of Distance 

Statistics. Geographical Analysis 1 992, 24:3. 

Stewart L, Gosling R, Griffin J, Gesase S, Campo J, Hashim R, Masika P, 
Mosha J, Bousema T, Shekalaghe S, Cook J, Corran P, Ghani A, Riley EM, 
Drakeley C: Rapid assessment of malaria transmission using age-specific 
sero-conversion rates. PLoS One 2009, 4:e6083. 
Stowers AW, Kennedy MC, Keegan BP, Saul A, Long CA, Miller LH: 
Vaccination of monkeys with recombinant Plasmodium falciparum apical 
membrane antigen 1 confers protection against blood-stage malaria. 
Infection and Immunit 2002, 70:6961-6967. 
Polley SD, Mwangi T, Kocken CH, Thomas AW, Dutta S, Lanar DE, 
Remarque E, Ross A, Williams TN, Mwambingu G, Lowe B, Conway DJ, 
Marsh K: Human antibodies to recombinant protein constructs of 
Plasmodium falciparum Apical Membrane Antigen 1 (AMA-1) and their 
associations with protection from malaria. Vaccine 2004, 23:718-728. 
Hardelid P, Williams D, Dezateux C, Cubitt WD, Peckham CS, Tookey PA, 
Cortina-Borja M: Agreement of rubella IgG antibody measured in serum 
and dried blood spots using two commercial enzyme-linked 
immunosorbent assays. Journal of Medical Virology 2008, 80:360-364. 
Gay NJ, Vyse AJ, Enquselassie F, Nigatu W, Nokes DJ: Improving sensitivity 
of oral fluid testing in IgG prevalence studies: application of mixture 
models to a rubella antibody survey. Epidemiology and Infection 2003, 
130:285-291. 

Vyse AJ, Gay NJ, Hesketh LM, Pebody R, Morgan-Capner P, Miller E: 
Interpreting serological surveys using mixture models: the 
seroepidemiology of measles, mumps and rubella in England and Wales 
at the beginning of the 21st century. Epidemiolgy and Infection 2006, 
1341303-1312. 

World Health Organization: Malaria programme reviews: a manual for 
reviewing the performance of malaria control and elimination 
programmes. World Health Organization;, Trial 2010. 
Mueller I, Galinski MR, Baird JK, Carlton JM, Kochar DK, Alonso PL, del 
Portillo HA: Key gaps in the knowledge of Plasmodium vivax, a neglected 
human malaria parasite. Lancet Infectious Diseases 2009, 9:555-566. 
Baird JK: Eliminating malaria - all of them. Lancet 2010, 376:1883-1885. 
Programme National de Lutte Contre le Paludisme: Directives nationales 
pour le diagnostic et le traitement du paludisme. Ministere de la Sante, 
Republique de Djibouti 2007. 

Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against 
malaria. Bulletin of the World Liealth Organizatio 2000, 78:1401-141 1. 
Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate 
and Plasmodium falciparum infection in African children. Nature 2005, 
438:492-495. 

Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic 
reproductive number for malaria and its implications for malaria control. 

PLoS Biology 2007, 5:e42. 

Smith DL, Drakeley C, Chiyaka C, Hay SI: A quantitative analysis of 
transmission efficiency versus intensity for malaria. Nature 
Communications 2010, 1:108. 

Health Metrics Network: Statistics save lives: strengthening country health 
information systems, draft. Geneva: Health Metrics Network; 2005. 
Gething PW, Noor AM, Gikandi PW, Ogara EA, Hay SI, Nixon MS, Snow RW, 
Atkinson PM: Improving imperfect data from health management 
information systems in Africa using space-time geostatistics. PLoS 
Medicine 2006, 3:e271. 

Cibulskis RE, Bell D, Christophel EM, Hii J, Delacollette C, Bakyaita N, 
Aregawi MW: Estimating trends in the burden of malaria at country 

level. American Journal of Tropical Medicine & hlygiene 2007, 77:133-137. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



64. 



Fox E, Bouloumie J, Olson JG, Tible D, LIuberas M, Shakib SO, Parra JP, 
Rodier G: Plasmodium falciparum voyage en train d'Ethiopie a Djibouti. 

Medicine Tropicale 1 991 , 5 1 :1 85-1 89. 

Carteron B, Morvan D, Rodhain F: The question of endemic malaria in 

Republic of Djibouti. Medicine Tropical 1978, 38:299-304. 

ICA meeting: Financing Transport for Growth in Africa - A Brief 

Memorandum of the Ethiopia-Djibouti Railway, [http://www.icafrica.org/ 

fileadmin/documentS/Transport_Meeting/S4-Djibouti-Ethiopia_Railway-Final- 

EN.pdf], Accessed December 2010. 

Federal Ministry of Health: The Ethiopia national malaria Indicator Survey 

2007. Ministry of Liealth, Federal Democratic Republic of Ethiopia 2008. 
Nyarango P, Gebremeskel T, Mebrahtu G, Mufunda J, Abdulmumini 0, 
Ogbamariam A, Kosia A, Gebremichael A, Gunawardena D, Ghebrat Y, 
Okbaldet Y: A steep decline of malaria morbidity and mortality trends in 
Eritrea between 2000 and 2004: the effect of combination of control 
methods. Malaria Journal 2006, 5:33. 

Snow RW, Alegana VA, Okiro EA, Gething PW, Patil P, Tatem AJ, Linard C, 
Sparrow A, Moloney G, Borle M, Yusuf FE, Amran J, Hay SI, Noor AM: 
Estimating the Plasmodium falciparum morbidity and mortality burden 
2005 and 2009 in Somalia: Combining models of population 
distribution, time-space changes in malaria infection risk and the 
epidemiology of malaria disease burden. Report prepared for UNICEF- 
Somalia 2010. 

Ministry of Public Health and Population: The Yemen national malaria 
indicator surveys of 2009. National Malaria Control Programme, Sana'a, 
Yemen 2010. 

Tatem AJ, Qiu Y, Smith DL, Sabot 0, All AS, Moonen B: The use of mobile 
phone data for the estimation of the travel patterns and imported 
Plasmodium falciparum rates among Zanzibar residents. Malaria Journal 

2009, 8:287. 

DEPHA (Data Exchange Platform for the Horn of Africa): [http://www.un- 
spider.org/sites/default/files/DEPHA.pdf], Accessed December 2010. 



Pre-publication history 

The pre-publication history for this paper can be accessed here: 
http://www.biomedcentral.eom/1471-2334/l 1/121/prepub 



doi:1 0.1 1 86/1 471 -2334-1 1-121 

Cite this article as: Noor et al:. Establishing the extent of malaria 
transmission and challenges facing pre-elimination in the Republic of 
Djibouti. BMC Infectious Diseases 201 1 11:121. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



